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ABSTRACT 


A Carousel  Remotely  Piloted  Vehicle  (RPV)  Recovery 
System  consists  of  a vertical  mast  which  supports  a horizon- 
tal member  that  is  free  to  rotate.  The  approaching  RPV  is 
caught  by  vertical  cables  suspended  from  the  horizontal 
member  and  the  kinetic  energy  of  the  RPV  is  dissipated 
through  the  motion  of  the  recovery  system.  This  thesis 
presents  a simplified  dynamic  analysis  to  describe  the 
motion  of  an  RPV  after  impact  with  a recovery  system  simu- 
lating the  carousel  system.  The  equations  of  motion  for 
the  RPV  were  obtained  from  Lagrange's  equations  and  modeled 
using  the  Continuous  System  Modeling  Program  on  the  IBM 
360-67  computer  at  the  V/.R.  Church  Computer  Center,  NPS. 

The  results  showed  that  damping  was  required  on  motion  of 
both  the  horizontal  member  ejid  the  suspended  cables  in  order 
to  prevent  possible  damage  to  the  150  poiuid  RPV  at  an  as- 
sumed 50  knot  approach  speed.  A sensitivity  analysis  was 
pei-formed  by  varying  the  system  design  parameters  and  the 
results  are  presented  in  the  body  of  the  thesis. 
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I.  INTRODUGTIOn 


Extensive  efforts  on  Remotely  Piloted  Vehicles  (RPV) 
have  been  underway  for  many  years,  covering  a v/ide  range  of 
applications  and  missions,  Ref.  1.  It  is  only  recently, 
however,  that  both  the  need  and  the  technological  capability 
for  mini-RPVs  have  coalesced  to  the  point  that  the  develop- 
ment of  these  systems  has  become  imperative.  A mini -RPV 
is  generally  defined  as  an  PlPV  that  v/eighs  less  than  200 
lbs  mission  loaded. 

An  essential  part  of  an  RPV  system  is  the  recovery 
system.  The  Army  has  recently  experienced  several  mini- 
RPV  losses  due  either  directly  or  indirectly  to  the  recovery 
method.  The  methods  used  to  recover  mini -RPV’ s vary  from 
catching  the  RPV  with  an  arresting  cable  similar  in  prin- 
ciple to  that  used  on  an  aircraft  carrier,  to  guiding  the 
mini-RPV  directly  into  a vertical  net  v/hich  is  designed  to 
flex  and  absorb  the  kinetic  energy  of  the  RPV.  A completely 
satisfactory  recovery  system  has  not  yet  been  established, 
hence  studies  into  alternate  recovery  systems  are  of  con- 
tinuing interest  to  the  Army. 

The  purpose  of  the  analyses  described  herein  was  to 
model  the  dynamics  of  one  alternate  recovery  system  candi- 
date, the  "carousel"  method.  The  carousel  method  was 
initially  described  as  a candidate  for  a recovery  procedure 
in  a concept  study  performed  by  Teledyne  Ryan  for  the  Army, 


11 


Ref.  2,  but  the  initial  concept  formulation  did  not  include 
significant  dynamic  considerations.  Therefore,  the  analyses 
conducted  in  this  thesis,  v/hich  are  primarily  oriented  to- 
wards system  dynamic  response  behavior,  should  serve  in  the 
role  of  providing  additional  information  for  concept  feasi- 
bility considerations. 

The  carousel  recovery  system,  shov/n  in  Figure  1,  was 
conceptualized  in  this  analysis  as  being  composed  of  a 
vertical  mast  that  supported  a horizontal  member  which  was 
allov/ed  to  rotate  about  the  vertical  mast.  Hanging  from 
the  horizontal  member  are  cables  that  will  be  caught  by 
the  v/ing  tips  of  the  approaching  RPV.  The  RPV's  kinetic 
energy  will  be  dissipated  by  the  recovery  system  through 
various  means  that  v/ill  be  evaluated  in  this  analysis. 

The  liPV  must  be  brought  to  rest  without  allowing 
e::cessive  RPV  swing  since  excessive  sv/inging  motion  coiild 
result  in  the  RPV  hitting  the  recovery  system  structure. 

This  consideration  v;as  of  paramount  concern  when  evaluating 
the  merits  of  recovery  system  parameter  variations. 

For  the  purpose  of  this  analysis  the  mini-RPV  to  be 
recovered  was  assumed  to  be  the  AQUILA,  which  is  a U.3. 

Army  mini-RPV  developed  by  Lockheed  Llissiles  and  Space 
Company,  Ref.  1,  and  is  sketched  in  Figure  2.  It  has  the 
following  nominal  characteristics:  wing  span,  12  ft; 
weight  150  lb;  approach  speed  50  knots;  maximxim  loading, 

6g  vertical  and  axial,  3g  lateral. 
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Figure  1 Carousel  Recovery  System 
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Figure  2 Aqulla  Mini-RPV 


II.  ANALYSIS 


A.  ASSUIviPTIONS 

Simplifyins  assumptions,  as  defined  below,  were  made 
when  establishing  the  equations  of  motion  using  the  pro- 
cedures suggested  by  Lagrange's  equations.  Ref.  3,  in  order 
to  establish  the  basic  principles  of  the  carousel  recovery 
system  dynamics.  Inclusion  of  system  details  such  as  re- 
covery boom  elasticity  and  RPV  mass  moments  of  inertia  were 
considered  as  items  which  might  tend  to  mask  the  basic 
question  of  feasibility  at  the  present  stage  of  analysis. 
Under  these  assumptions,  the  results  supported  the  viability 
of  the  cg,rousel  recovery  method  and  a basis  was  established 
upon  which  further  analyses  may  be  performed  when  the  sim- 
plifying assumptions  are  removed. 

The  recovery  system  was  initially  assximed  to  be  con- 
servative, that  is,  friction  losses  v/ere  ignored  and  as 
soon  as  the  RPV  contacted  the  recovery  system  the  engine 
\7as  secured.  This  assumption  was  subsequently  changed  by 
the  addition  of  damping  to  the  recovery  system  configuration. 
The  horizontal  member  was  assumed  to  be  rigid.  A careful 
choice  of  material  and  structural  characteristics  for  the 
horizontal  member  can  minimize  its  dynamics  and  allow  it  to 
be  approximated  by  a rigid  member.  The  cables  hanging  from 
the  horizontal  member  were  repla.ced  by  massless  thin  re- 
covery rods.  The  RPV  '.vas  treated  as  a point  mass  and  was 
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captured  by  one  recovery  rod,  Treating  the  RPV  as  a point 
mass  simplified  the  coordinate  system  used  in  the  analysis 
by  allowing  the  mass  moment  of  inertia  of  the  RPV  about  its 
principal  axes  to  be  neglected.  The  RPV  was  also  assumed  to 
approach  perpendicular  to  the  plane  formed  by  the  horizon- 
tal member  and  the  vertical  recovery  rod  that  caught  it. 

The  approach  velocity  of  the  RPV  was  assumed  to  transfer 
undiminished  to  the  recovery  rod  giving  the  rod  an  initial 
angular  velocity  only  in  the  0 direction.  Reference  should 
be  made  to  Figure  3 where  the  orientation  of  the  generalized 
coordinates  e,  0 and  ij/  is  portrayed.  These  coordinates 
are  a form  of  Euler  angles,  independent  of  each  other,  that 
together  completely  specify  the  position  of  the  RPV  relative 
to  the  fixed  reference  frame  x-y-z.  First  0 defines  the 
rotation  of  the  moving  coordinate  system  x'-y'-z'  about  the 
fixed  axis  z.  It  describes  the  rotation  of  the  horizontal 
member  in  the  x-y  plane.  Next  0 defines  the  motion  of  the 
RPV  about  X*  within  the  moving  coordinate  system.  This  is 
portrayed  in  Figure  3 as  the  motion  of  the  RPV  projected  on 
the  y'-z'  plane.  Finally,  ij/  depicts  the  remaining  motion 
of  the  RPV  after  it  rotates  through  0 and  0 . 

The  system  identified  in  Figure  3 is  the  reference  or 
base  configuration  that  was  assumed  for  the  analysis  based 
on  judgement  of  the  RPV's  weight  and  size.  In  this  refer- 
ence configuration  the  horizontal  member,  R^»  vvas  28  ft 
long  and  made  of  6-inch  a.luminxim  I-beam  having  a mass 

p 

moment  of  inertia  of  1159  ft-lb-sec  about  pivot  "0", 
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The  recovery  rod  R^,  was  6 ft  long.  The  RPV  of  mass  4.66 
slugs  (150  I’os)  was  suspended  at  the  lower  end  of  R^.  »Vith 
an  assumed  approach  speed  of  50  knots  (84.4  ft/sec),  the 
initial  angular  velocity,  0^,  of  the  massless  rod  R2  v/as 
14  rad/sec  according  to  the  relation  of  equation  (1): 

1 Approach  velocity  (ft/sec) 

= a'  (hji  ' ^ (1) 

A linear,  torsional  spring,  of  spring  constant  K,  lo- 
cated at  "0",  the  pivot  of  the  horizontal  member,  was 
initially  used  to  resist  the  rotation  of  the  horizontal 
member,  i/ithout  any  damping  in  this  reference  configura- 
tion it  was  further  assumed  that  the  spring's  stored  energy 
could  be  eliminated  from  the  system  through  a ratchet  re- 
lease mechanism,  after  the  horizontal  member  had  swung 
through  its  initial  maximum  displacement.  The  ratchet 
mechanism  was  not  modeled,  as  it  was  felt,  and  later  con- 
firmed by  the  results,  that  the  critical  motion  of  the  RPV 
would  occur  shortly  after  impact  with  the  recovery  system. 
Simulation  of  the  ratchet  mechanism  would  have  no  effect 
on  these  results. 

B.  OVERVIEW  OP  AlULYSIS 

The  general  approach  of  the  analysis  was  to  model  the 
RPV  motion  after  engagement  with  the  recovery  system  using 
the  generalized  coordinate  system  0-0- indicated  in 
Pig.  3.  If  the  assumptions  of  a rigid  horizontal  member 
and  a point  mass  RPV  representation  were  removed,  a minimxim 
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of  five  more  generalized  coordinates  would  have  been  re- 
quired. The  equations  of  motion  for  the  three  degree  of 
freedom  system  were  obtained  using  Lagrange's  equations 
for  a conservative  system.  ?rom  these  equations  the 
acceleration,  velocity  and  position  of  the  RPV  as  a function 
of  time  were  obtained  in  component  form  along  each  coordi- 
nate direction  through  the  use  of  the  Continuous  System 
I/Iodeling  Program  (CSIfP).  The  input  parameters  were  then 
varied  and  their  effect  on  resulting  dynamic  motion  studied. 

Changes  to  the  basic  recovery  system  which  included  the 
addition  of  damping  terms  were  later  entertained  as  a con- 
sequence of  initial  results. 


C.  DETAILED  ANALYSIS 

1.  Lagrange  * s Equation 

The  equations  of  motion  were  established  using 
Lagrange's  equation,  Ref.  3.  Three  equations  of  motion 
were  obtained  for  the  system,  since  there  were  by  assump- 
tion only  three  generalized  coordinates:  the  coordinates 
e,  0 and  , as  defined  in  Fig.  3.  Lagrange's  equation 
in  the  conservative  form,  without  any  generalized  forces,  is: 


i ( 

« Jc.i 


i—  1»  2,  3> 


(2) 


= generalized  coordinate  (qQ_=  0 »q2=0  >CL2=l/r) 
q^  = time  derivative  of  q^  coordinate 
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T = T(qj_,43_)  = Kinetic  energ5’-  (ft-lb) 
V = V (q^)  = Potential  energy  (ft-lb) 


The  kinetic  energy  was  defined  relative  to  the  inertial 
reference  frame  by  the  following  expression: 


T 


6 + a 


(3) 


where 


I 


0 


i.I 

e 


^RPV 


rioment  of  Inertia  of  the  horizontal  member 
about  pivot  point  "0"  (ft-lb-sec^) 

Mass  of  RPV  (slugs) 

Angular  velocity  in  0 coordinate  direction 
(rad/sec) 

Velocity  expression  for  RPV  derived  in  the 
next  section  (ft/sec) 


The  potential  energy  included  the  gravitational  term  acting 
upon  the  RPV  mass  and  the  stored  energy  in  the  torsional 
constraint  spring.  Equation  4 expresses  the  potential 
energy  relative  to  the  inertial  frame  of  reference  as 
follows: 

V = MgR2(^l-cos0cosT^yi-K6^  (4) 

where 

R2  = length  of  recovery  rod  corresponding  to  RPV 

distance  from  the  horizontal  member  following 
recovery,  (ft) 

K = spring  constant  of  spring  at  pivot  of  horizontal 
member  (ft-lb) 


20 


Since  its  exactness  is  critical  to  the  evaluation  of 


Lagrange’s  equations  of  motion,  the  RPV  velocity  is  treated 
separately  in  the  ensuing  section  using  the  generalized 
coordinates. 

2.  Derivation  of  RPV  Velocity 

To  derive  an  expression  for  the  RPV's  velocity 
after  impact  with  the  recovery  system,  the  position  of  the 
RPV  within  the  moving  x'-y'-z’  coordinate  system  as  a func- 
tion of  the  generalised  coordinates  0,  was  transformed 
to  the  fixed  inertial  reference  system  x-y-z.  The  time 
derivative  of  the  RPV's  position,  expressed  as  a function 
of  the  three  generalized  coordinates  0,0  and  and  ref- 
erenced to  the  fixed  inertial  system,  represented  the  ab- 
solut-e  velocity  of  the  RPV. 

The  RPV's  position  within  the  moving  frame  x'-y'-z' 
v;as  exoressed  as  follows: 


'rt 


f _ 


sin  l/r 

R2  sin0cosT/^ 
-R2  COS0 cos  y/- 


(5a) 

(5b) 

(5c) 


Transformation  to  the  fixed  inertial  coordinate 


stern 

nay 

z'  ■>. 

X 

» = 

y 

COS0  - sin0 
sin0  cos  0 


r 


-R^sin"^ 

R^slnCpcos^l/' 


r 


+ < 


ij^cosQ 


z = 


= z'  = -R^cosCp  coslfr 


R-,sin0 

. 1 J (6a) 


(6b) 
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In.  an  alternate  form,  the  three  components  expressing  the 
absolute  position  of  the  RPV  mass  following  recovery  en- 
gagement may  be  expressed  as: 

X = R^cos0 -RgSin^cos  0 -R2Sin0  cos'j//'sin0  (7a) 

y = R^cos0  -R2sinl//^sin0 +R2sin0  cos  "^^cos  0 (7b) 

z = -R^cos0cos  ij/'  (7c) 

The  absolute  velocity  of  the  RPV  was  obtained  by  taking  the 
first  time  derivative  of  equation  7,  thereby  yielding: 

X = ^-R^sin0 +R2( sin‘^sin0 -sin0cos'^cos  e)]  e 

-(R2COS0  cos  'l/rsin0  )0  +R2(  sin0  sinl^sin  Q - 
cosT/Tcos  0 ) (8a) 

.y  = ^R^cos0 -R2(sin0"cos0 +sin0  cos’^sin0  ^ 0 

+ (R^cos0cos'l/rcos 0 )(p  -R2^cos0'sin0 +sin0sinl/r 
COS0  (8b) 

r • * — 

s = Rp  |(sin0cos  cos^sinl^)  ij/'j  (8c) 

The  expression  desired  for  use  in  Lagrange's 
2 

equation,  Vj^p^  , was  obtained  by  squaring  the  components 
separately  then  adding  them;  i.e., 

Vppv^  = x^+y^+s^  (9) 

-2-2  ’2 

The  expressions  for  x , y and  z are  equations  10a,  b 
and  c . 
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= ^-R^sin0  +R^(  sinl/z^sin  Q -sin0  cosl^cos  e)]"e= 

+R2^cos  0cos^'^sin‘' 0 0 

+R2^(sin0sinl/rsin0 -cos’0'cos0 
-2R2Cos0cos1^sin0  |-R^sin  0+R2(sinY^sin0  -sin0 

cosT/rcos0  ^00 

+ 2R2(  sin0Gin'0'sin0  -cocl/^cos  Q ) j^-R^  sin0  h-R^C  sinyr 

• • • 

sin0 -3in0cosl/rcos0  ) 0 ilf 

O • • 

-2E2‘=coe0  coslj/'slnQ  (sin0  slnlj/'slnQ  -coalj/'cosQ  )0'^ 

(10a) 


jR^cos0  -R2( sirxl/rcos0  +sin0 cosT^sin 0 )J  ^ 0^+R2^ 

2 2 2*2 
cos  0COS  1^  cos  0 0 

+R2^  jcosT^sin  0 +sin0sin'0'cos  0j 

+ 2R2COS  0cosl^cos  0 |r^cos0  -R2(sin0'cos  0+sin0 
cosl/Osin0  )J0  0 

-2R2  jconYrsin0 +sin0cin'^//'cos  0j  j^R^cos0 -R2(  sinY/" 
COS0  +sin0cos Y^sin 0j  Qij/' 

-2R2^co30cos  Y^'ccs  0 |cos  Y^sin0  +3±n0s±n^ 
cos  0I  1^0 

-I  ^ ^ (10b) 
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• P PrP  2 2 2 2*2 

pin  0 cos  "^0  +003  0 sin  Y/" l//^ +2sin0cos0 

sinl/^cosl//"  01^'']  (lOc) 

•2  -2  *2 

'.Vhen  X , y , z were  combined,  several  cancellations 

2 

occurred  and  the  final  expression  for  Vpp.^  simplified  to: 


^RPV^  = |^R^^+R2^  sin^1^+R2^sin^0  cos^'^-2RpR2Sin'l/rj 

+R2^oos'^  Tp- 0^+R^^  ■\j/^+2R^cos0  cosijriR^-R^slnllr)  Q 0 
+ i-2R2_^2^^^0  sinl/r+2R2^sin0  ) 0 ij/' 

(11) 

3 . Derivation  of  Equations  of  Llotion 
a.  0 Equation  of  Motion 

2 

Now  that  an  expression  for  V^py  has  been 
obtained,  Lagrange's  eo^uation  can  be  evaluated  according  to 
section  II.G.l.  The  equation  was  evaluated  for  each  gener- 
alized coordinate  e , 0 and  Ip'  . The  0 equation  of  motion 
v/as  obtained  by  evaluating  the  Lagrange's  equation  using 
the  0 generalized  coordinate. 


d /^T  \ 21  + - 0 

^ 30  ' je  36"  ° 


Again; 


Kinetic  energy,  T-I-Iq  0 ^+i-IvIVj^PY^ 

Potential  energy,  V=MgR2  [l-cos0  cos  4-jK  0 ^ 


21  _ 
30 


0 


(12) 


(13) 

(14) 

(15) 
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(16) 


^^2|j=i^0  +M  sln^ 0 cos^ljr -ZR^R^slnljr)  e 

+R2Cos0cos'|/^(R]_--R2  sin  1/^)0  +(-R]_R2sin0  sim0+ 
R2^sin0)  I/'' 

2 * 2 2 * 

+ (2R2  sin'0'coc'0'i^+2R2  sin0cos0cos  ij/' 0 

-2R2^sin^0  cosT/rsin10 T^-2R^R2Cos-|0' t^)  0 

« 2 

-R2sin0cos'l/r(R^-R2Sin'0')  0 -R^slnlfrcos  0{R^-R^ 


sinl/^)  0 

2 . '•  ■ . ' ' 

-R2  cosT,^'cos0 cos'0''l|^0  -R1R2CP  s0  sin-0'0  ■0'-RiR2 

• 2 

sin0cos'0''0' 

+R2^cos0  01^  (17) 


Substituting  equations  15 » I6  and  17  into  equation  12  and 
dividing  through  by  (1.1^2^)  yields  the  following  ©equation 
of  motion; 
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I R 

— ^+(-|T=-sinl/r)^+cos^l/rsin^0  ) 0 

R, 

+sin.(^(l  ^—sinT^)  ij/' 

+ COS0  (^cos'|^-^sin2l/r) 

2 ■ 

+ sin2l//' cos  0 -2^cosT/^)l/re 


2 * * 

+ sin2<pcoG  Ip'  0 Q 


"1 

+COG0  (-2^sinT/r-cos2^+l)  (p 
R 

1 • 2 

-Gin  COG  ^-:>sin2l/r)  0 


R 


;-2 


- — Gin0coG  Yn/r' 


^ r r.> 


2=0 


(18) 


b.  0 Equation  of  I, lot  ion 

Similarly  the  0 equation  of  notion  v/aG  ob- 
tained from  the  following  Lagrange 'g  equation: 


0. 0 


d f 3T  ' 

15>, 


^ = 0 

^0  ()0 


(19) 


whe  re 
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R2^cos^l/r0  +R2COG0  cosl/rCR^-RgSinl/r)  0 


Q • • • • 

-2R2^co  sl/rsirL"^  1^0 -R2Sin0  cos  R^-RgSin  l/r)  0 0 

-R2Sim^co30  (R^-R2Sinl//‘)0i^  -R2^cos0  cos2^T^0j 

(20) 


= r.i[R2^sin0cos0cos^-^0^-R2Sin0cosl/r(R^-R2Sini^)00 


2 * * 

+ (-R^R2Cos0  sinT/''+R2  cos0)0'^j 


^ = r.IgR«sin0cosl/r 
90  2 y 


(21) 


(22) 


■ Substituting  equa,tions  20,  21  and  22  into  equation  19  and 

2 

dividing  through  by  (ll?2  ) yields  the  follov/ing  0 equation 
of  motion; 


Rl 

7^cos0  coslf/'--hcos0  sin21^ 
'^2 

+ 003“!/^  0 


-sin2'l/^  ijrcp 
-COG0  ( COS2'^+l)  0 

-^s±n20  cos^lfr  0^ 


+ ^sin0cosl/h=  0 


e 


(23) 
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c.  Equation  of  Motion 

The  Ip"  equation  of  motion  was  similarly  obtained 
using  the  following  Lagrange's  equation: 


d(  2l\  - + 


(24) 


wne  re 


d / ^T 


at 


2 r V ^1  * ' ^1 

= M[T/4+(-  g-cos0sim^+cos0  )00  - ^cosl/rsin0-^0 


" "I 

+ (-  ^sin0sinT^+sin0  )0J 


(25) 


R. 


R^^M  [( sinl/rcos'l^-sin^0  cos'\p‘s±n\f/~-  •^cos'i^) 

* 2 T_  2 

-cos'^sim/r0  +(-  ^cos0s±n'\J/{R^-R^s±n'\l/')-CQS  ij/- 


COS0)  0 0 


R. 


_ ^sin0cos'^  0 


S:Cos0sinY^ 


2I 


(26) 


(27) 


Substituting  equations  25,  26  and  27  into  equation  24  and 

2 

dividing  through  by  the  (I.ER2  ) .yields  the  following  l/r 
equation  of  motion: 
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R 

(-  ^sin0sim/r+sin0  )0  + ijr 

+ COS0  (l+cos2y^)  Q 0 

+ (^sin2‘^(r(  sin^0 -1)+  ■^cos'l^)  0^ 

2 

• O Q. 

+o"sin2'l/^  (2)  + #-cos(2)  sinT//'=  0 

^ ^2  ^ ^ (28) 

4.  Solving  the  Equations  of  IJotion 

The  equations  of  motion  were  solved  for  time  his- 
tory information  of  the  RPV's  position,  velocity  and  accel- 
eration using  the  IBM  Continuous  System  Modeling  Program 
(CoMP)  which  v/as  available  as  resident  software  on  the  IBM 
360-67  computer  located  in  the  V/,R.  Church  Computer  Center, 
Naval  Postgraduate  School.  The  use  of  CSr.IP  as  a high-level 
programming  tool  is  quite  similar  to  the  familiar  engineer- 
ing language  of  FORTRAN,  and  programming  details  may  be  found 
in  the  test  by  Speckhart  and  Green,  Ref.  5.  Details  of  the 
actual  program  used  are  presented  in  the  Appendix.  To  use 
this  software  package,  the  three  equations  of  motion  were 
first  solved  for  the  highest  order  derivative  0 , 0 and  ’if' 
in  terms  of  the  three  generalized  coordinates  and  their  first 
time  derivatives.  Then,  in  a manner  similar  to  an  analog 
computer  solution  of  differential  equations,  the  CS?;EF  program 

proceeded  to  integrate  these  equations  numerically  for  values 

for  6,0,1/^,  0,0  and  T/T,  which  were  fed  back  into  the 

••  ••  • • 

equations  for  0 , 0 and  ’ij/'  . This  iterative  method  of  , 
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solution  was  repeated  at  each  time  increment.  The  resulting 
values  5or  0 , 0 and  were  listed  as  a function  of  time 
and  from  these  the  motion  the  RPV  was  studied  to  assess  the 
feasibility  of  the  recovery  system  when  a particular  set 
of  structural  parameters  and  initial  RPV  conditions  v/as 
specified. 

Since  the  governing  equations  of  motion  (eqns.  18, 

23  and  2S)  were  inertially  cross  coupled,  it  was  necessary 


to  rearrange  the  equations  into  an  inertially  uncoupled  form 
that  v.'ould  be  acceptable  for  solvition  by  the  CSI'.IP.  In  matrix 
form,  equations  18,  23  and  28  may  be  expressed  as; 


ail  a^2  ^13- 

CD: 

f ^1  h 

^21  ° ^23 

< 

t 

> = < 

b2 

^■31  ^32  ^ 

0 > 

■w  b^  ^ 

v.'he  re 


^11 

ai2 


*13 


-23 


I-  S p'  2 2 

^ -sin'^/r)  +COS  -^sin  0 


2 ' ^ ^ 2 


a^i  = sin0(l isin"^) 


R. 


^21  " cos  0(^003*^ -j-cin  T/r) 


cos^  "[jr 


^32  “ 


1 


(29) 

(30a) 

(30b) 

(30c) 

(30d) 

(30e) 


30 


(30f) 


= ^sin0cos'^'^+sin0  (^cosY/"-:Vsin2  ■^)  (p‘^ 


-(sin  cos"^0  -2^cosT^)0  -sin20  0 0 

-.(l-cos2l/''-2Tpsin'^)  cos  0-^0-  t? 


= ^coG‘'T/rsin2  0 0^+cos0  (l+cos2l/^)  00  +sin2^'I^0- 


^sin0cos^ 


(30-) 


= ({;sin2 ■l/'^cos^0  - ■^cos'l^)  0^-;Vsin2'|//'0^  -cos0(l+ 


cos2l/^)00  - ^cos0sin'|^ 


(30h) 


Sub 3 1:1  tut ing  equations  30a  through  3 Oh  in  matri:-:  expression 
29»  0>  0 - and  ijjr  were  solved  using  Cramers'  n.ile,  Ref.  4. 

0 =1:^1  T^=l'^2l  0 = l:^l 

I A|  (31a)  ^ I A|  (310)  I A]  (31c) 


vvnere 


ail 

p 

^'12 

P 

..13 

®-21 

0 

^'23 

= ^12  ^23  ^31"^11  ^23'*'^21  ^13 

^31 

^32 

0 

(32a) 

°1 

ai2 

^13 

^2 

0 

^23 

= ^3  2-12  ^23“^1  ^23'^^2  ^13 

1 

0 

(32b) 

31 


- 


i'H  = 


^11 

^13 

^21 

^23 

- ^3i(^3_^23“^2^13^”^'^. ^^11^23” 

^21®-13^ 

^31 

^3 

0 

(32c) 

^11 

^12 

aoi 

0 

b2 

= 2-2]_(b2aQ^2)“(623-3_]_“D-i^a2]_)+b2 

®-31 

1 

’“3 

( ■"^21^12  ^ 

(32d) 

2quations  iOa  ‘through  32d  were  progranuned  on  the  IBIJ  360-67 
computer  using  CSI.IP,  the  details  of  v/hich  are  described  in 
the  Appendix.  The  prograra  evaluated  0»0,^r,  0,0,1^, 
e,  0 and  ijr  S.S  a.  function  of  tine  from  the  initial  con- 
dition of  RPV  inpact.  ‘The  evaluation  of  the  computer 
generated  time  history  results  is  described  in  the  next 
sectio'n. 
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III.  SYSTEIJ  DESIGIT  iUIAlYSIS  AI^ID  RESULTS 

A.  UI\T[)AI.IPED  SYSTELi 

The  reference  or  base  configuration  of  the  recovery 
system  described  in  section  II. A.  and  Fig.  3 was  the  first 
configuration  considered  in  the  analysis.  From  this  con- 
figuration each  parameter  was  varied  individually  for  various 
flPV  approach  speeds  to  assess  the  system' s sensitivity  to 
their  changes. 

From  the  outset  of  the  investigation  it  was  felt,  as 
supported  by  early  results,  that  the  motion  of  the  RPV 
along  the  <p  coordinate  would  be  a primary  system- constraint . 
A guideline  value  for  0maximum  ( 0max)  . was -chosen  to  be  90 
degrees,  v/hich  corresponds  to  the  recovery  rod  (R2)  reaching 
the  horizontal.  For  the  recovery  system  to  pass  its  first 
requirement,  an  RPV  approaching  at  50  knots  v/ould  have  to 
be  recovered  v/ith  the  maximum  amplitude  of  0 not  exceeding 
90  degrees. 

The  results  of  modeling  the  reference  configuration 
suggested  that  energy  disipation  would  have  to  be  included 
to  bring  0ma::  do'/.n  to  90  degrees  or  less.  V/ithout  any 
damping  the  system  exhibited  unacceptably  high  magnitudes 
for  0 that  corresponded  to  the  RPV  completely  rotating 
about  the  horizontal  member  v;hen  approach  speeds  exceeded 
25  knots.  This  can  be  seen  by  referring  to  Table  I where 
0 is  listed  for  the  base  configuration  as  a function  of 
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VARIATION  FOR  UNDAMPED  SYSTEM 
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tine  after  inpact  for  approa,cii  speeds  ranging  iron  7 knots 
(00=2  rad/sec  Ref.  Ea_uation  1)  to  50  knots  ( 00  =14  rad/sec). 
Rote  that  for  0q="  rad/sec,  or  approximately  a 25-lcpot 
approach  speed,  the  magnitude  of  (p  at  0.7  seconds  after 
impact  v.'as  greater  than  130  degrees.  The  RPV  has  started 
to  rotate  over  the  horizontal  member.  Variation  of  the 
system  parameters  did  not  correct  this  deficiency. 

3.  DA-,IFIKG  ADOED  TO  THE  HORIZONTAL  iIELSER 

In  an  effort  to  correct  the  excessive  R?V  sv/ing  in  the  0 

coordinate  direction,  a Rayleigh  dissipative  energj'  term, 

3,  v.’as  introduced  into  Lagrange's  equation: 

^ + ll  . 12-0 

- . (33) 

• n 

v/here  3=vG-,  @ 

This  additional  term  precipitated  through  the  ecj_uation 

of  motion  derivation  of  section  II. C.  v/ith  the  addition  of 

> 

only  a 0^6  term  to  the  0 equation  of  notion,  equation  18. 
This  change  v/as  easily  incorporated  into  the  computer  pro- 
.gram  as  indicated  in  the  Appendix. 

All  initial  estimate  of  the  magnitude  of  the  damping 
constant  nas  obtained  using  the  simplified  linearised 
pendulum  equation  below.  The  damping  value  was  varied  to 
detect  the  system's  sensitivity  to  G^. 
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C C2;.= 


id  0.15, 

v/hich.  v;c.c  selected  to  introduce  nust  a sr.all  amount  of 


for  values  of  1^=1159  Id-ft-sec  , IC=50  fo-ld 


damping,  tiie  initial  estimate  of  ‘‘^^.s  75  ft-lb-sec. 

Tlie  results  of  adding  damping  to  the  pivot  of  the 
horizontal  member  in  the  base  confi^ration  indicated  that 

it  v/as  only  partially  successful  in  re  due  ins  0mam.  i'his 
can  be  seen  by  referring  to  fable  II  where  0 is  listed  for 
this  configuration  as  a function  of  time  after  impact  for 
approach  speeds  ranging  from  7 Icnots,  v/here  0q=2  rad/sec, 
to  50  knots,  where  0^=14-  rad/sec.  bote  that  the  magnitude 
of  0 v/as  still  approaching  180  degrees  for  approach  speed 
initial  conditions  in  emcess  of  35  knots  ( 0^=10  rad/sec). 

Variation  of  system  pa.rameters  including  did  little 
to  further  reduce  0max,  as  shown  in  Figures  4 through  8, 

v/here  the  base  configuration  of  11^=28  ft,  22=6  ft,  .',1=4.66 
slugs,  1^=1159  ft-lbs-sec^,  K=50  ft-lb  and  0^=75  ft-l'o-sec 
is  plotted  with  only  the  indicated  parameter  varied,  ceteris 
pa.ribus . 

G.  OX.IPIKG  ALQdb  TO  HKGO’/FRY  ROD 

In  a.  further  effort  to  reduce  the  maximu-m  swing,  0ma;c, 
of  the  recovery  system,  damping  was  a.lso  added  to  the 
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Figure  4 Effects  of  Damping  on  0[nax 
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Figure  7 Effects  of  Horizontal  Member  Length  R 
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recovery  rod.  The  added  damping  was  included  in  the  system's 
model,  without  cross  coupling,  oy  modifying  the  Hayleigh 
dissipative  energy  function  to  the  form  of: 

(35) 

The  additional  dissipation  term  precipitated  through  the 
equation  of  motion  derivation  described  in  section  II. C.  as 
the  addition  of  a 0^0  term  to  the  0 equation  of  motion, 
equation  23.  This  change  was  also  easily  incorporated  into 
the  computer  program  as  indicated  in  the  Appendix. 

An  estimate  of  C2  was  obtained  in  a manner  similar  to 
that  described  in  the  previous  section  using  the  simple 
pendulum  equation; 

0 +2  ^ (0^0  + (jQ^0=  0 

v/here 


Only  a slight  amount  of  damping,  02=3.24  ft-lb-sec, 
added  to  the  recovery  rod  of  the  base  configuration,  in 
addition  to  the  horizontal  member  damping,  was  sufficient 

to  reduce  0max  to  approximately  90  degrees  at  RPV  approach 
speeds  of  50  knots  ( 0o“^^  rad/sec).  This  can  be  seen  by 
referring  to  Table  III  where  0 is  listed  as  a function  of 
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time  after  impact  for  approach  speeds  ranging  from  7 knots 
( 00=2  rad/sec)  to  50  loiots  ( 0^=14  rad/sec).  The  maxirainn 

magnitude  of  0 and  'ij/'  for  the  modified  conf iguration  at  the 
50-knot  RPV  approach  speed  v/as  50  degrees  and  63  degrees, 
respectively. 

Further  study  of  the  modified  base  configuration  was 
performed  to  obtain  the  magnitude  of  the  RPV*s  linear 
acceleration.  The  absolute  value  of  the  RFV*s  linear 
acceleration  v/as  obtained  using  the  following  equation; 

/ 2 » 2 

Acceleration  = y x +y  +z  (37) 

where  x,  y and  s are  defined  by  equations  7a,  b and  c.  A 
plot  of  equation  37  as  a fiinction  of  time  for  the  modified 
base  configuration  is  presented  in  Figure  9.  It  ma^’’  be  ■ 
noted  that  the  magnitude  reduces  to  below  one  "g"  within 
one  second  after  the  RFV  contexts  the  recovery  system. 

The  sensitivity  of  0max  to  each  of  the  system  para- 
meters is  indicated  in  Figures  10  through  16,  where  the 
base  configuration  of  R^=28  ft,  R2=6  ft,  i.i=4.66  slugs, 
Iq=1159  ft-lb=sec2,  k=50  ft-lb,  C^=75  ft-lb-sec  and  0^= 

3.24  ft-lb-sec  is  plotted  with  only  the  indicated  parameter 
varied,  ceteris  paribus.  It  can  be  noted  from  Figure  12 
that  the  system  was  sensitive  to  the  recovery  rod  length, 

R2»  which  may  be  interpreted  as  an  implication  that  the 
engagement  point  for  RPV  capture  on  the  vertical  cables 
has  considerable  importance  to  system  dynamic  response. 
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Figure  10  Effects  of  Danq>ing  Cj^  on  0n,ax  with  C2  Damping 


2.00  ft- lb 
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Figure  11  Effects  of  D 
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Figure  13  Effects  of  RPV  Mass 
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Figure  14  Effects  of  Horizontal  Member 
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Figure  15  Effects  of  Moment  of  Inertia 


The  systeni  was  also  sensitive  to  the  value  of  C^,  but 
nore  sensitive  to  C2.  This,  however,  should  not  present  a 
problem,  as  and  vvill  be  set  by  the  design  and  will 
not  vary  in  operation. 

The  system  was  relatively  insensitive  to  variations  of 
the  point  along  the  horizontal  member,  where  the  RPV 

was  caught.  The  system  was  also  relatively  insensitive  to 

« 

variations  in  the  moment  of  inertia,  of  the  horizontal 

member  and  the  mass,  !VI,  of  the  RPV.  The  maximum  displace- 
ment of  0 v/as  not  affected  at  all  by  changes  in  the  spring 
constant,  K. 

D.  i:,IPL2r.ISNTING  DALiPING  OR  A RECOVERY  CABLE 

The  clamping  desired  at  the  pivot  of  the  horizontal 
member  can  be  readily  implemented  in  actual  construction. 
However,  it  was  not  readily  apparent  how  the  damping  of 
the  vertical  rod,  which  was  a cable  in  the  actual  system, 
could  be  implemented. 

Two  simple  pendulum  designs  were  analysed  and  modeled 
using  the  same  techniOj.ues  described  to  study  the  recovery 
system.  The  intent  v/as  to  see  if  pivotal  damping  could 
be  applied  effectively  to  control  the  swing  of  a cable. 

The  first  design  is  shov/n  in  Figure  17  and  simulates 
a cable  design  for  the  vertical  massless  recovery  rod  that 

was  used  in  the  recovery  system  analysis.  V/ith  no  pivotal 
damping,  an  initial  velocity  was  applied  to  II2  "to  malce 

approach  a high  angle  without  exceeding  I80  degrees. 
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Then  damping  v/as  applied  to  the  pivot  of  the  cable  support 
r.I^  to  judge  its  effect  on  reducing  the  magnitude  of  (p^. 

The  results  are  plotted  in  Figure  18.  This  approach  of 
applying  pivotal  damping  did  not  effectively  reduce  p2m'^x 
as  desired. 

Another  design,  shov/n  in  Figure  19,  v/as  also  analysed 
to  see  if  it  could,  through  pivotal  damping,  control  the 
maximum  magnitude  of  0^.  Again  the  design  simulated  a 
cable  instead  of  the  massless  recovery  rod  used  in  the 
analysis.  A spring  v/as  used  at  the  pivot  to  give  a 

near  horizontal  static  deflection.  With  no  damping  applied, 
an  initial  velocity  was  found  for  M2  that  caused  02m3.x 
approach  I80  degrees.  Pivotal  damping  v/as  then  applied 
.to  to  observe  its  effect  on  the  magnitude  of  02* 
results  shov/n  in  Figure  20  indicate  that  pivotal  damping 
applied  to  this  pendulum  design  v/as  also  ineffective  in 
reduo  in-;  0^^.,. 

./hile  neither  of  these  designs  demonstrated  effective 
control  of  0^^  they  by  no  means  exhaust  the  possible 
methods  of  implementing  damping  in  a cable  design,  and 
further  study  is  possible. 
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IB.OB 

Damping  Constant  (ft- lb-sec) 

Figure  18  Results  of  Pendulum  A Analysis 


IV.  CONCLUSIONS  AifD  RECO?4I.IENDATIONS 


For  th.e  carousel  recovery  system  to  be  feasible  at 
approach,  speeds  of  50  knots,  a relatively  small  amoimt  of 
damping  v/as  required  at  the  pivot  of  the  recovery  rod  in 
addition  to  damping  at  the  pivot  of  the  horizontal  member. 
The  damping  at  both  pivots  v/as  required  to  prevent  the 
maximum  swing  of  the  point  mass  RPV  from  exceeding  approx- 
imately 90  degrees  in  the  0 direction,  the  coordinate  that 
exhibited  the  largest  change  in  magnitude. 

The  recovery  system  was  sensitive  to  the  vertical 
height  of  the  RPV  at  the  time  it  was  captured.  The  greater 
the  distance,  R^,  between  the  RPV  and  the  horizontal  member, 
the  greater  the  swing  of  the  RPV  relative  to  the  horizontal 
member.  V/hen  selecting  parameter  values  such  as  the  damp- 
ing constants  under  the  most  severe  conditions,  the  RPV 
should  be  assumed  to  have  contacted  the  recovery  system  at 
the  maximiim  length  of  the  recovery  rod  or  cables. 

The  recovery  system  was  relatively  insensitive  to 
variation  in  the  distance  the  RPV  was  from  the  vertical 
mast  when  it  was  caught.  The  system  was  also  relatively 
insensitive  to  variations  in  the  mass  of  the  RPV  and  the 
moment  of  inertia  of  the  horizontal  member. 

Variations  in  the  spring  constant  for  the  torsional 
spring  located  at  the  pivot  of  the  horizontal  member  had 
little  effect  on  the  swing  of  the  RPV.  Since  it  was  not 
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important  to  controlling  the  motion  of  the  RPV,  it  could  be 
removed  from  the  design. 

The  magnitude  of  linear  acceleration  exhibited  by  the 
RPV,  with  the  recovery  system  design  that  incorporated  both 
horizontal  member  and  recovery  rod  damping,  was  within 
Aquila  load  limits. 

Based  on  this  analysis  the  carousel  recovery  system  is 
still  viable;  however,  it  is  apparent  that  the  system  will 
be  sensitive  to  approach  techniques  used  during  the  recovery. 

It  is  recommended  that  this  analysis  be  continued,  in 
order  to  treat  the  RPV  as  a distributed  mass  with  its  ovm 
moment  of  inertia  and  coordinate  system,  and  caught  by  tv/o 
recovery  rods  or  cables,  one  for  each  wing  tip,  instead  of 
one.  Expan(?ing  the  scope  of  the  recovery  system  model  would 
provide  a Lire  detailed  analysis  of  the  RPV's  linear  lateral, 
longitudinal  and  normal  accelerations  and  hence  insure  that 
they  were  individually  within  the  RPV’s  design  limits. 

It  is  also  recommended  that  the  effects  of  axi  RPV 
approaching  at  angles  to  the  recovery  system  be  investi- 
gated; this  should  result  in  increased  motion  in  the  'ij/' 
direction. 

It  is  further  recommended  that  continued  study  should 
be  made  of  a means  of  applying  damping  to  the  swing  of  the 
vertical  recovery  cables. 


APPEIOIX 


:iODSLIITG  THE  EQUATIONS  OF  MOTION  USING  CSr.IP 

An  application-oriented  software  package  called  the 
Continuous  System  Modeling  Program  (CSMP)  was  used  to  model 
and  solve  the  nonlinear,  cross-coupled  equations  of  motion. 
GSI.IP  is  a powerful  software  package  that  can  handle  this 
type  of  problem  routinely.  In  the  analysis  performed  in  this 
thesis  two  integrations  using  a fourth  order  variable  inte- 
gration step  Hunge-Kutta  technique  on  each  of  the  three 
equations  of  motion,  Ref.  5.  This  sophisticated  integration 
technique  has  the  advantage  of  automatically  adjusting  the 
time  increment  of  integration  to  meet  the  demands  of  the 
dynamic  conditions  of  the  simulation.  The  absolute  value  of 
the  estimated  integration  error  and  the  relative  magnitude 
of  the  estimated  error  are  compared  with  user  specified 
error-bounds.  The  step  size  is  then  adjusted  to  meet  the 
desired  error  criteria.  In  the  analysis  performed  in  this 
thesis  the  default  value  of  .0001  for  both  the  absolute  and 
relative  errors  was  used. 

The  program  structure  of  GSMP  is  composed  of  three 
segments;  Initial,  Dynamic,  and  Terminal.  Generally,  data 

statements  will  appear  in  the  Initial  segment.  In  addition, 
calculations  that  are  required  to  be  performed  only  one 
time  during  the  simulation  can  be  conveniently  placed  in 
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this  segment.  The  Dynamic  segment  is  usua,lly  composed  of 
structural  statements  that  describe  the  dynamics  of  the 
system  or  explicitly  describe  a set  of  differential  equa- 
tions. The  Terminal  segment  is  the  last  segment  in  the 
program  and  is  usually  made  up  of  output  control  statements. 

The  actual  program  used  in  this  analysis  for  the  \xn- 
daraped  system  will  now  be  described.  List  of  symbols  used 
in  the  Computer  Programs; 

R1  - Length  of  horizontal  member  in  feet 

112  - Length  of  recovery  rod  in  feet 

M - Ivlass  of  RPV  in  slugs 

I.IIO  - IJoment  of  inertia  of  horizontal  member  in  ft-lb- 
sec^  about  its  pivot  point 

K - Spring  constant,  in  ft-lb/rad,  of  spring  at  pivot 
of  horizontal  member 

X5I  - Initial  angular  velocity  of  the  recovery  rod  in 
the  0 direction  in  units  of  rad/sec  also  equal 
to  approach  velocity  of  RPV  in  ft/sec  divided 
by  R2 

p 

G - Acceleration  of  Gravity  32.2  ft/sec 

XI  - e 

X2  - 0 

X3  - T/d 

X4  - e 

X5  - 0 

x6  - ijr 

X4D0T  - S 
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X5D0T  - 0 

X6D0T  - T/r 

CO  - G^,  damping  constant  in  ft-lb-sec  for  the  damper 
located  at  the  pivot  of  the  horizontal  member 
C02  - Cgj  damping  constant  in  ft-lb-sec  for  the  damper 

located  where  the  recovery  rod  attaches  to  the 
horizontal  member 

AD  - Determinant  of  A as  defined  in  Section  II. C. 4 

ADI  - Determinant  of  A1  as  defined  in  Section  II. C. 4 

AD2  - Determinant  of  A2  as  defined  in  Section  II. C. 4 

AD3  - Determinant  of  A3  as  defined  in  Section  II. G. 4 

B1  - is  b^  defined  in  Section  II. C. 4 

B2  - is  b2  defined  in  Section  II. C. 4 

B3  - is  b^  defined  in  Section  II. G. 4 

ill  - is  a^2_  defined  in  Section  II. G. 4 

A2  - is  defined  in  Section  II. C. 4 

A3  - is  defined  in  Section  II. C. 4 

A4  - is  3.^2^  defined  in  Section  II. G. 4 

A,  B,  C,  D,  E,  F,  H,  1,  S,  T - were  used  to  simplify 

program  notation  and  are 
defined  in  the  program 

A description  of  the  program  shovm  in  Figure  21,  follows; 
Liiies  1-3:  flake  up  the  initial  segment,  parameters 

that  would  remain  constant  during  the 
integration  were  assigned  here. 

Lines  4-34:  Make  up  the  Dynamic  segment. 

Lines  5-13:  Simplifying  notation. 
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Figure  21  Computer  Program 


Lines  14-28:  The  algebra  used  to  solve  for  0,  0 and 'ij/' 
prior  to  integration  as  described  in 
Section  II. G. 4. 

Lines  29-31:  Perfox-m  the  first  integration  from  angular 
accelerations  to  angular  velocities. 

Lines  32-34:  Perform  the  second  integration  from  angular 
velocities  to  angular  position. 

Lines  35-41:  Ivlahe  up  the  Terminal  segment. 

Line  36:  Limits  the  time  period  that  will  be  modeled. 

Lines  38-39:  Prints  output  as  a function  of  time. 

V/hen  damping  was  added  to  the  horizontal  member  the 
only  changes  that  occurred  to  the  program  v/ere  the  follov/ing: 


Line 

2: 

add 

,C0  = 

75 

Line 

23:' 

add 

1 

o 

o 

* 

X4 

V/hen  damping  v/as  also  added  to  the  recovery  rod  the 
additional  changes  that  occurred  to  the  program  were  as 
follows: 

Line  2:  add  ,002  = 3«24 
Line  24:  add  -C02  * X3 
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